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We show that supersymmetry and R-symmetry can be dynamically broken in a long- 
lived metastable vacuum of SQCD with massive and massless flavors. The vacuum re- 
sults from a competition of a (leading) two-loop effect and small "Planck" suppressed 
higher-dimension operators. This mechanism provides a particularly simple realization of 
dynamical SUSY and R-symmetry breaking, and as such it is a good starting point for 
building phenomenologically viable models of gauge mediation. We take a preliminary step 
in this direction by constructing a complete model of minimal gauge mediation. Here we 
find that the parameters of the model are surprisingly constrained by the hidden sector. 
Similar mechanisms for creating long-lived states operate in a large class of models. 
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1. Introduction 



Dynamical supersymmetry breaking (DSB) via strong gauge dynamics has long been 
an attractive mechanism for explaining the large hierarchy between the weak scale and 
the Planck scale [Q. Early, seminal models of DSB were important "existence proofs" - 
however, they tended to be rather complicated (see e.g. 0-[5| for reviews and references). 
This was tied to the assumption that the DSB vacuum had to be the global minimum of the 
potential. By relaxing this assumption and allowing for metastability, much simpler models 
of DSB have recently been constructed, starting with the discovery of metastable vacua 
in SQCD with massive flavors ||. (Of course, metastability and its connection to SUSY 
breaking are not new ideas; for earlier works see e.g. [[7|-|TT|.) This work has opened new 
model building avenues, especially for models of direct gauge mediation. Since metastable 
SUSY-breaking models can be vector-like, they can naturally have large flavor symmetries. 
This makes them more amenable for direct gauge mediation, where a subgroup of the flavor 



symmetry is gauged and identified with the SM gauge group |T^| . 

Unfortunately, despite all its positive features, the metastable vacuum of SQCD with 
massive flavors has one major phenomenological defect: it preserves a continuous U(1)r 
symmetry. Without breaking this symmetry somehow, the gauginos of the visible sector 
would remain massless.0 

Various approaches have been tried to get around this problem. Motivated by 
the fact that the low-energy effective theory of massive SQCD is simply a generalized 
O'Raifeartaigh model, more sophisticated O'Raifeartaigh models with R-charge assign- 
ments different from or 2 were studied in jjTJ], and it was shown that such models could 
lead to spontaneous R-symmetry breaking. A detailed study of the phenomenology of 



such models was undertaken in [[1511 ; however, no UV complete dynamical realization of 
these models is known. (Non-renormalizable realizations based on massive SQCD with 
baryonic deformations were constructed and studied in |T6|JT7[1 .) Alternatively, one can 
gauge an auxiliary £7(1) symmetry of the problem; the competition between the gauge 
and Yukawa couplings can lead to a spontaneous R-symmetry breaking minimum of the 



1 Two comments. First, the R-symmetry is broken explicitly by nonperturbative effects but 
this is generally too small to allow for realistic gaugino masses. Second, we are assuming here 
standard Majorana gaugino masses. An interesting alternative whose phenomenology has recently 



been revisited in [13] is to have Dirac gaugino masses, in which case the R-symmetry need not be 
broken at all. 
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one-loop effective potential [Tc-20]. However, here there generally has to be some fine 
tuning of the couplings in order to achieve the desired R-symmetry breaking ]21[] . Finally, 
a third (and perhaps most obvious) possibility is to break the approximate R-symmetry 
explicitly, but on general grounds |22[ this generically leads to proximal SUSY vacua and 
some tuning has to be invoked (see e.g. J23|-|25| for models). One possibility that has been 
explored is to have explicit R-symmetry breaking only through Planck-suppressed cou- 
plings to messengers [^6,^7|. However, in this case the messenger scale and the messenger 



SUSY splittings come from different sources, so there are in general extra CP phases in 
the gaugino masses when doublet /triplet splitting of the messengers is taken into account. 
Similarly, messenger-parity is not guaranteed and there is a danger of inducing spontaneous 



breaking of the SM gauge group [28 



In this work we present another mechanism of spontaneous R-symmetry breaking 
which can be viewed as an improvement upon previous attempts, inasmuch as it is simple, 
more natural, dynamical and phenomenologically viable. The crucial observation is that 
all the abovementioned attempts at spontaneous R-symmetry breaking assumed that the 
R-charged pseudo-moduli obtain their potential at one-loop. In some instances, however, 
pseudo-moduli do not obtain their potentials until two-loops, and then the story can be 
quite different. Here, we will study one of the simplest illustrations of this phenomenon: 
SQCD with both massive and massless flavors. In this model, the meson composed of the 
massless quarks remains massless at one- loop in the low-energy effective theory |29| . Much 
above the strong coupling scale the theory exhibits nonperturbative runaway behavior; 
below the strong coupling scale, the two-loop potential was recently computed in the low- 
energy effective theory and was shown to be a monotonically decreasing function of the 
meson VEV [30|. Thus the model probably has runaway for all values of the field and has 
no vacuum at all. 

In this paper we show how the runaway is automatically stabilized by "Planck" sup- 
pressed higher-dimension operators in the superpotential. The result is a local vacuum far 
from the origin but well below the cutoff scale of the effective low-energy theory, where 
both SUSY and R-symmetry are spontaneously brokenJl 



A similar idea in a different context was proposed in [31| to explain the hierarchy between 
the Planck scale and the GUT scale. See also [^9,32,33] for recent models based on massive 
SQCD where dangerous pseudo-moduli directions were stabilized with higher-dimension operators. 
However in these latter examples the pseudo-moduli in question were stabilized at the origin and 
were not the sources of R-symmetry breaking. 
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Aside from the SUSY-breaking minimum, the higher-dimension operators also restore 
a supersymmetric solution. We show that the distance between the supersymmetric and 
SUSY-breaking solution is parametrically large, ensuring that the metastable vacuum is 
parametrically long-lived. 

Consequently, SQCD with massive and massless flavors has all the necessary ingre- 
dients, including large spontaneous R-symmetry breaking, for SUSY model building. We 
illustrate this by writing down the simplest application of our model, namely a complete 
model of "minimal" gauge mediation |8|,[M|]. The role of the singlet coupling to messen- 
gers in our model is played by the meson composed of the massless electric quarks. As 
we couple the meson to the messengers through non-renormalizable interactions 



m 



in the UV. We find that if the runaway is stabilized by a quartic coupling in the UV, 
surprisingly stringent constraints on such models of minimal gauge mediation come from 
the requirements of calculability in the SUSY-breaking sector together with non-tachyonic 
messengersJl As we show further this model can be successfully retrofitted |36| and all the 



couplings can be explained in the context of discrete gauge symmetries |37]. 

Of course, the real promise of our new mechanism of dynamical SUSY and R-symmetry 
breaking is in its potential applications to models of direct gauge mediation, where there 
is no invariant distinction between messenger and SUSY-breaking fields. Since SQCD 
with massive and massless flavors is vector-like, it automatically comes equipped with 
large global flavor symmetries. It is easy to imagine gauging a subgroup of these flavor 
symmetries and identifying it with the SM gauge group. Then SUSY breaking will be 
automatically communicated to the SM via gauge mediation, and unlike in , the MSSM 
gauginos will acquire viable soft masses. In this context one still has to overcome the 
long-standing challenge of the Landau pole problem, since there will generally be lots of 
extra matter charged under the SM gauge group. This is a very interesting problem that 
we will reserve for future work. 

Our paper is organized as follows. In section 2 we review SQCD with massive and 
massless flavors. We pay special attention to the theory far out along the pseudo-moduli 
space (where the runaway will ultimately be stabilized) . In this regime, the theory can be 
described by an effective supersymmetric Lagrangian, and the leading- log approximation to 



3 These constraints can be relaxed if one chooses to stabilize the runaway by even higher- 
dimension operators. The main drawback of this possibility is that a more complicated discrete 
symmetry is needed to maintain such a structure. 



3 



the runaway potential has been calculated as part of a more general analysis in [38], using 



the techniques of wavefunction renormalization [ p9| , |40| . However, to keep the discussion 
in this paper self-contained we present a sketch of the calculation in the appendix. In 
section 3 we introduce the stabilization of the two-loop runaway potential via higher- 
dimension operators, discuss the spectrum of particles around the SUSY-breaking vacuum, 
and derive constraints on the parameters of the model which follow from calculability and 
longevity. Section 4 is dedicated to constructing and analyzing a model of minimal gauge 
mediation using massive+massless SQCD as the SUSY-breaking sector. 



2. SQCD with Massive and Massless Flavors 

2.1. An overview of the model 



In this section we review the main results of ||30|| , while simultaneously establishing our 
notation and conventions (which differ slightly from those of ||30|1 ). We will consider SUSY 
QCD with SU(N C ) gauge group and Nf flavors in the free magnetic phase (N c < Nf < 
3N c /2). We will take Nfo < N c of the flavors to be massless with the other Nfi = Nf —Nfo 
having equal mass.@ When all the flavors are massive (but light compared to the strong 
coupling scale A), this model exhibits a metastable supersymmetry breaking state @. 
However, when some of the flavors are massless, the situation is very different: some of 
the pseudo-moduli of the low-energy effective theory are not stabilized at one-loop f2S 



Instead, they acquire a potential at two-loops and it was recently shown in [|3(| that this 
potential is monotonically decreasing along the pseudo-moduli space, leading to runaway 
behavior. 

Now let us describe the low-energy effective theory in more detail. Using Seiberg 
duality |[41|| , we have a weakly coupled description at low energies in terms of an IR-free 
SU(N) gauge theory (with N = Nf — N c ). The matter content consists of a gauge singlet 
Nf x Nf meson matrix $ and Nf flavors of magnetic quarks (p, (p. It will be convenient 
in the following to split the Nf flavor indices into Nfo and Nfi sized blocks, i.e. 

-ft *::)• *-(£)• *--(£)• <"> 



4 It is trivial to ensure the naturalness of this structure by imposing appropriate (discrete) 
symmetries which allow mass terms for some of the quarks but not others. We will have more to 
say about this in the next sections. 
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with an Nfi x Nfj matrix and fi, iff Nfi x N matrices. The fields can be normal- 
ized to have canonical Kahler potential plus uncalculable higher-dimensional corrections 
suppressed by powers of A. In this normalization, the superpotential of the theory is 



W = hTr^ijifjifi - /i^ 2 Tr$n 



(2.2) 



where h is an 0(1) coupling and the traces are taken over the uncontr acted flavor indices 
of the different fields. In what follows we will assume without loss of generality that all 
the couplings are real and positive. Note that although the global symmetry is no longer 
SU(Nf), we are still denoting all the Yukawa couplings by the same symbol h - at the 
scale A this is the case up to corrections which vanish as £ —>■ 0. Of course, below the scale 
A these couplings will run differently, and we will take this into account when necessary. 

The representations of the fields under all the gauge and global symmetries preserved 
by the superpotential (|2.2|) are listed in the table below. Note that there are four indepen- 
dent U(l) symmetries (including the R-symmetry), neglecting anomalies as in 0. 





SU(N) 


SU(N f0 ) L 


SU(N f0 ) R 


SU(N fl ) 


U(1)b 


U(l)i 




U(1) F 
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1 


1 
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□ 
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1 
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□ 


1 











<Po 
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1 
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-1 


1/2 


<h 
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1 


□ 


-1 











<Po 


□ 


1 


□ 


1 


-1 


-1 





1/2 



(2.3) 

The reason for the funny £7(1) charges will become apparent in the next section. 

Since Nfi = Nf — Nfo > Nf — N c , rank conditions mean that SUSY is sponta- 
neously broken at tree- level, as in §. The tree- level scalar potential is minimized along 
the following pseudo-moduli space: 



$i 







X 
$01 = (0 Y) , 

with vacuum energy 



$io= (° 
arbitrary $oo 



fo = ifo = 

2, ,4 



<Pl 



xx = 



N 



V = (N fl - N)h 2 » 



(2.4) 
(2.5) 



5 



Here x, \ are N x N matrices, X is an (Nfi — N) x (Nfi — N) matrix, Y and Y T 
are (Nfi — N) x Nfg matrices, and $00 is an Nfo x Afyo matrix. At the origin of the 
pseudo-moduli space ( |2.4| ), the entire gauge symmetry and part of the global symmetry 
are broken: 

SU(N) x \sU(N fl ) x U(1) B ] -> [s?7(7V)£, x SJ7(iV)i - AT) x £7(1)' B 1 , (2.6) 

with the other global symmetries remaining unaffected. 

As mentioned at the beginning of this section, the main novelty introduced by the 
massless quarks is the fact that at one-loop, only some of pseudo-moduli are lifted by 
the Coleman- Weinberg potential p9fl . Specifically, all the fields except $00 are lifted and 
stabilized at the origin by the one-loop potential. So to understand the dynamics of this 
model we are driven to a two-loop computation on the pseudo-moduli space parameterized 
by $oo- This has been considered in detail in |3(J. There it was found that the origin 



of field space is destabilized, in other words, $00 acquires a negative mass squared term 
around the origin 

V eff = V (^) 2 N(N C - N f0 ) (l + y - log 4^ Tr($J $oo) + , (2.7) 



, 2 

where ah = j~- Moreover, it was shown in [30 that the two- loop potential for $00 is 



monotonically decreasing, leading to runaway behavior. This result has important physical 
implications for ISS based model building and it can be interpreted geometrically in the 



appropriate brane configuration [30]. In the next subsection we will study the model in 



more detail in the ($00) H regime. 

2.2. Effective theory at large ($00) 

Since the pseudo-modulus $00 has a runaway potential, we are motivated to study the 
theory in the large $00 ^ A* regime. There the analysis of the pseudo-modulus potential 
simplifies considerably - we can use the wavefunction renormalization techniques of [j3~9| , |40|| 
to compute the two-loop potential in the leading-log approximation using only one-loop 
anomalous dimensions. 

To be more precise, we are interested in the theory for generic ($00) i n the regime 



hp < h($ao) < A . (2.1 
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The point is that below the scale h($>oo) we can integrate out supersymmetrically the Nfo 
magnetic quarks coupling to $00 • Then we are left with an effective supersymmetric theory 
of the form 

(2.9) 

W eff = hTr$u(pi(pi - h^Tr^n - hTr ¥>i$io$oo $oi¥>i , 
where the wavefunction factors in the Kahler potential are functions of $00 and where 
. . . includes both the uncalculable 0(1/ A) corrections from the duality together with the 
calculable C(l/($oo)) corrections from integrating out (po, ipo. 

The effective theory still has an SU(N) gauge symmetry; the gauge coupling can 
be either IR or asymptotically free depending on whether Nfi > 3N or Nfi < 3N, 
respectively. The latter would be problematic because then the Kahler potential would no 
longer be well-behaved around the origin where we need to compute loops to understand 
the SUSY-breaking vacuum. Thus in order to have a reliable SUSY-breaking vacuum we 
must require 

N f i > 3N = 3(N f - N c ) . (2.10) 

This goes beyond the usual free-magnetic phase requirement Nf > 3N. Note that ( |2.10| ) 
is equivalent to 

iV/o < 3N C - 2N f < N c . (2.11) 

The effective theory is essentially the "macroscopic model" of || but with a wave- 
function factor in the Kahler potential that depends on $00 • The theory still breaks SUSY 
via the rank condition and stabilizes all the pseudo-moduli. The vacuum energy is: 

V Q = (N fl - N)h 2 ii A Z^l . (2.12) 

This can be interpreted as the effective potential for $ Cb after using the RGEs to compute 
Zq,^ in the leading-log approximation. The computation of Z<$, xx has been performed 
as a special case of a more general survey of pseudo-moduli. Here we simply quote 
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the answer; we will reproduce the calculation in appendix A for the sake of completeness. 

lo g ^iv(^) 2 Tv( log ?^)\..., (2.13) 

where . . . are subleading in the two-loop leading-log approximation. Plugging this 
into ( |2.12| ), we have in the leading-log approximation 

V e ff = -(N fl -N)Nh^ (^) 2 Tr flog^j • (2.14) 

Of course, (|2.14|) can also be computed more directly by taking the large field limit of the 
two-loop potential computed in [[30] . 
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3. Stabilizing the Runaway 

3.1. Adding a small quartic term 

In this section, we will show how the runaway can be turned into a metastable SUSY- 
breaking vacuum by allowing for generic non-renormalizable superpotential interactions. 
One can easily imagine that such terms are present at the Planck scale, generated by the 
theory of quantum gravity. Alternatively, they could be trivially generated at an interme- 
diate scale by integrating out massive matter in some renormalizable UV completion. In 
any event, the point is that, since the runaway direction is a dimension two meson in the 
electric theory, the lowest-dimension non-renormalizable operators - dimension four inter- 
actions amongst the electric quarks - become mass terms for the runaway fields in the IR. 
So the runaway is automatically stabilized by generic higher-dimensional superpotential 
interactions. 

Instead of considering the most general quartic terms for the electric quarks, we will 
focus only those involving the Nfo massless quarks; these become in the IR: 

SW = ^rn figif2 g 2 (^ 00 )f igi (^oo)f 2 g 2 . (3.1) 

It would be interesting to study in detail more general deformations involving $oi 5 $io 
and $n, but we will not do so in this paper. We will briefly discuss the possible effects 
of such deformations at the end of the next subsection. The upshot is that, while such 
deformations may complicate the details of our analysis, it is not obvious that they will 
change the qualitative conclusions. 

For Nfo > 2, the quartic deformations in ( |3.1| ) necessarily break the SU(Nfo)L X 
SU(Nfo)n symmetry down to a subgroup. For simplicity, we will choose to preserve 
SU(Nfo) o- We will further simplify the analysis by considering only one of the two 
independent deformations that respect this symmetry: 

SW = -/le/xTr $q . (3.2) 
2 



One could of course consider an arbitrary mixture of (3J2) and the other SU(Nfo)o sym- 
metric deformation (Tr$oo) 2 - However, since this complication will have no effect on the 
essential physics that we will discuss below, we will ignore it. Notice that in ( |3.2| ) we have 
written the deformation parameter in units of /i/x, with a new dimensionless parameter e. 
This will simplify the notation throughout. 
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In order for the deformation to not overwhelm the negative two-loop mass-squared 
for $00 ( |2.7|) around the origin, we need e to be a small parameter (as its name suggests), 
specifically 

« * £ ■ (3-3) 
Ordinarily this small parameter would seem rather artificial; however, here the UV com- 
pletion can make it natural - all we need is that the UV scale that suppresses ( |3.2|) in the 
electric theory is sufficiently large compared to A. 

We are interested in whether the deformation can stabilize the runaway at large $00 ^ 
(x. In this regime, the effective theory (|2.9|) deformed by (|3.1|) is a valid description. Thus 
the vacuum energy along the pseudo-moduli space $00 is 

V = (Nfi - N)h 2 ^Z^ + /i 2 eV 2 Tr $ 00 $ o ■ (3.4) 

Here we are neglecting the wavefunction renormalization of $00 itself - since we are going 
to be balancing the two terms in ( |3.4p against each other to find the metastable vacuum, 
any corrections from Z$, 00 are obviously going to be subleading. 

Substituting the calculation (|2.14| ) of the previous section, it is straightforward to 
minimize the potential, 



Veff = -(Nfi - N)Nh 2 ^ (^) Tr ( log ) + $00*00 , (3.5) 

and find the local minimum at 



oo) « £ yiog i l Nf0 , - 5 - yJ{N fl - N)N (g) \ . (3.6) 

The corrections to Q are of the form ($ o> -> ($00) x (l + O L^, los £° S 5 S) )) and 
they come not only from minimizing the potential (|3.5|) , but also from the subleading- log 
corrections to the wavefunction Z$ tl . It is very important (and fortunate) for our purposes 
that with the leading-log potential alone, we can capture the location of the $00 VEV ( |3~^| ) 
up to small corrections. 

In the regime e <C the minimum of the potential is indeed located at ($00) /U, 
which ensures that the approximation used to obtain the minimum is self-consistent.il 



5 Note that this is a stronger assumption than ( ^.3[) . In fact, ( |3.3|) is sufficient for the existence 
of a minimum away from the origin - since the potential turns up at large fields, as long as the 
origin is destabilized such a minimum must exist. In this paper we will only consider the more 
restricted case e <C ^ where the effective theory is valid. 
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Another requirement comes from demanding $qo ^ A, so that the effective description in 



terms of the magnetic theory can be trusted. According to (|3.6|) , this implies 



47T A 

(In fact, we will see below that an even stronger bound is required for calculability.) Since 
these requirements (|3.3| ), ( |3.7|) can obviously be satisfied for the proper choice of param- 
eters, we conclude that a viable minimum of the potential can be obtained by balancing 
the two- loop runaway potential against a tree- level mass term for $00 • 

3.2. Properties of the metastable vacuum 

Now let us describe the properties of the metastable vacuum in more detail. The 
deformation (|3.2|) explicitly breaks some of the global symmetries ( [2.3|) . The remaining 
symmetries are: 

SU(N) x \sU(Nf ) D x SU(N fl ) x U(1) B x Z7(l)i x U(1) R ] . (3.8) 



Note that the choice of U(l) charges in (|2.3j ) makes the symmetry breaking here especially 
simple. 

In the metastable vacuum, the symmetries are broken spontaneously by ($00) 7^ 
and (pi, cpi 7^ to 

SU(N) x \sU(N f0 ) D x SU(N fl ) x U(1) B x f7(l)i x U(1) R 

(3.9) 

SU(N) D x SU(N f0 ) D x SU(N fl - N) x U(1)' B x U(l) 1 

Thus we see the crucial difference with ||: the U(1)r symmetry is completely broken in 
the metastable vacuum! This can have nice applications to phenomenology, as we will see 
in the next section. 

There is a rich spectrum of fluctuations around the metastable vacuum. We will now 
describe these, roughly in order from heaviest to lightest. 

1. First of all, the quark superfields (po, Lp$ are very heavy with C(/i($oo)) masses. At 
the tree-level, the squarks have small SUSY splittings dm 2 ~ /i 2 e/x($oo)- 



2. From the description in terms of the effective theory (|2.9| ), ( |3.4| ), we see that the 
spectrum below the scale h($>oo) basically factorizes into an ISS-like spectrum for 
$11, (pi, (px phis the spectrum of the additional modes $01, ^10, $oo- The fields in 
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the ISS sector have masses ~ h\i (except for the pseudo-moduli, which have one-loop 
suppressed masses- squared, and the Goldstone bosons). Some of the fields have large 
dm 2 ~ h 2 fi 2 splittings. Since these fields have been given a detailed description in ||, 
we will move on to describing the new fields. 

3. Let us split up $01 and $10 into two submatrices each: 



where A, A T are N x Nfo and Y, Y T are (iV/i — N) x Nfo. Then from the third term 
in the effective superpotential (|2.9|) , we see that A, A acquire masses ~ /i / u 2 /($oo) 
from the VEV of (p v They also have SUSY splittings 5m 2 ~ ^F$ 00 ~ In 
order for these to not make the scalar components of A, A tachyonic, ($00) cannot 
be too large: ($00) At A- Fortunately, this is guaranteed by the formula for the 
met ast able vacuum (|3.6| )! 

4. Meanwhile Y, Y are massless at tree-level. Since they have a suppressed coupling 
~ h[i/(<fr 00 ) to the SUSY-split "messenger" fields in cpi, <pi, their one-loop mass- 
squared will be similarly suppressed, m 2 ~ ( ) 2 (h[i) 2 . An explicit calculation 
shows this mass-squared is positive. Regarding fermions in this sector, the mass term 
ipY^y * s anowe d by the unbroken symmetries and it is indeed generated at one-loop. 
As before, to generate this operator we have to go through propagators which are 
suppressed by l/($oo)- 

5. Finally, $00 itself obviously has a mass ~ hen, except for Im(Tr$oo) which is the 
massless Goldstone boson of U(1)r breaking.! 

Finally, let us comment on the potential effects of other quartic operators besides the 
ones considered in ( |3.1| ). From the discussion of the spectrum, it becomes clear that adding 
operators of the type e/i$n$io, e//$n$oi to the superpotential will change our vacuum 
significantly. The reason is that the F-term of $n gives rise to a tadpole for the light Y 

6 Actually, the effective scalar theory ( |3.4| ) appears to possess the full SU(Nfo) 2 symmetry 
even though the superpotential ( |3.2| ) breaks it down to SU(Nfo)D- So when ($00) 7^ it would 
seem there should be many more tree-level massless modes than just the Goldstone boson of 
U(1)r breaking. In fact, these accidental massless modes are all lifted by considering the most 
general 5C/(A r /o)_D-symmetric deformation, instead of just (|3.2|) . 





(3.10) 
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modes, and shifts them far away from the origin. These operators might not change the 
qualitative conclusions, but for our analysis to strictly apply they have to be absent for 
some reason, e.g. by a SU(Nf Q )o x SU(Nfi) global symmetry or discrete symmetries like 
those we consider in section 4. By the same reasoning, operators of the form e/x$n$oo are 
also dangerous, so they should also be forbidden by symmetries. In fact, these operators 
are automatically forbidden by any symmetry which allows a linear term in $n, while 
simultaneously forbidding $00 and allowing $q . 

3.3. Calculability, nonperturbative effects and lifetime 

Finally, let us briefly address various issues related to the consistency of our analysis 
of the metastable vacuum: calculability, irrelevancy of the nonperturbative effects, and the 
lifetime. 

First, because the tree-level mass of $00 is so much smaller than /i, we need to be 
careful that uncalculable corrections to the Kahler potential do not overwhelm the tree- 
level superpotential for $oo- In particular, terms like 

in the Kahler potential, with c an uncalculable 0(1) number, contribute to the mass of 
$00 via F$ u ~ hfi 2 : 

Veff 3 ^(V 2 ) 2 $ 00 $oo • (3.12) 
In order for this to be negligible compared with the tree- level mass (Zie/u) 2 , we clearly need 

6 » £ . (3.13) 

Note that this is a stronger requirement than (|3.7|) which we deduced above by demanding 
that ($00) A. We will get back to this constraint in more detail in section 4. 

Next, the issue of whether nonperturbative effects can be neglected around the 



metastable vacuum boils down to the IR-freedom requirement (|2.10|) , just as in ||. In 
more detail, when $00 an d $11 have generic VEVs, then by the standard scale matching 
argument the dynamical superpotential is 

l/N 
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The calculability requirement (|2 . 1 ) also guarantees that the non-perturbative superpo- 



tential is an irrelevant (dim> 3) operator around $n = even for $oo 7^ 0. Thus it can 
be safely ignored when computing loops around the non-SUSY vacuum. 

Finally, let us address the issue of the lifetime. For this we need to know where the 
SUSY vacua are. At this point we will take h = 1 to simplify the expressions, because 
the value of h does not affect our analysis of the lifetime. Taking the ansatz $oo oc Ijv /0 
and $n oc IjVfj (and all other fields zero), it is straightforward to set the F-terms of the 
magnetic superpotential to zero (including (|3.14| )) and solve. Introducing the additional 
small parameter 

?=^-, (3.15) 



this yields (ignoring phases) 



3N c -2N f +N c -N fQ Nfi+N 



($00)SUSY =€ 2N ^ N fO e 2iV c -iV /0 ^ = ~2N c -N f0e 2N c -N foA 

2(3N C -2N f) 2(3JV c -2JV^)-JVj Q 
/ rf\ \ ~ 2iV c -iV fn , 2iV c -JV fn ,, 



(3.16) 



So we always have 

/u < ($oo)susy < A , ($u)susy > /" • (3.17) 
Since (&n)non-suSY = 0, (|3.17|) implies that the SUSY vacuum is always separated 



by an amount A$n ^> fi from the non-SUSY vacuum. This is enough to ensure a para- 
metrically long-lived metastable state, provided that the barrier also scales with A$n. 
The following argument shows that this must be the case. As in 0, the most efficient 
path between the non-SUSY and SUSY vacua is to first climb up to a local extremum with 
q = q = and V = Nfi^. Then the potential is schematically (now setting A = 1) 



V ~ N fl 



N f0 /N^N fl /N-l a 



2 



^00 Q ll - V + iV /0 ^00 Q ll _ e V® 



-N f0 /N-l^N fl /N 



'00 



2 



(3.18) 



Now the question is how far does one have to go in field space before the potential decreases 
to V non -suSY = (Nfi — N)fx 4 7 Since the vacuum energy comes primarily from the first 
term in ( |3.18| ), at the very least we have to decrease it by making 



Consider now all values of $oo an d $11 which satisfy ( |3.19| ). We claim that of all such 



values, either $ o >/ior <& n > \x. Indeed, if $u < /i, then $ o ^ n( 3N ~ N f^/ N fo and one 
can check that this satisfies $oo 3> {&oo)non-suSY , (&oo) susy ^> fi- Therefore the width 
of the potential barrier is guaranteed to be parametrically larger than //, either through 
A$n > /j or through A$qo ^ ^- This in turn ensures a parametrically long lifetime for 
the metastable vacuum. 
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3.4- Summary of the requirements on the metastable vacuum 

Since the discussion has been technical at times, we would like to conclude this sec- 
tion by collecting all the different requirements that need to be satisfied in order for the 
metastable vacuum to exist. 

In addition to requiring that the theory is in the free magnetic phase (N c < Nf < 
|iV c ), we also must require 

iV /0 < 3N C - 2N f Nfi > 3(N f - JV C ) (3.20) 

in order for the non-perturbative effects to be negligible around the metastable vacuum, 
and to ensure a parametrically long lifetime. 

The parameter e = m^> 00 /hfi must satisfy the constraints: 

f « e< ^ . (3.21) 

The lower bound comes from demanding that uncalculable Kahler potential corrections 
can be neglected relative to the tree-level mass term. It also ensures that ($oo) A. The 
upper bound comes from requiring that the tree-level mass term does not overwhelm the 
two-loop potential and stabilize the metastable vacuum at the origin. We will be assuming 
a stronger bound e C ^ throughout the paper so that the leading-log approximation to 
the effective potential ( |3.5| ) can be used. 



4. A Simple Model Building Application 

4-1. The model 

While the main emphasis of our paper is to present a new mechanism of metastable 
supersymmetry-plus-R-symmetry breaking, in this section we would also like to briefly 
touch upon how this mechanism can be applied to building viable and natural models of 
gauge mediation. (A more detailed study of the model building applications, in particular 
to models of direct gauge mediation, will be reserved for a future publication.) To that end, 
we will present here the simplest application of our mechanism to model building: using 
our SUSY-breaking sector in a model of "ordinary" or "minimal" gauge mediation |]8|,^4|. 
Although this may seem like a trivial application, we will see that various consistency con- 
ditions and phenomenological requirements force the model to be surprisingly predictive. 
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Our SUSY-breaking sector consists of the theory we studied in great detail in the 
previous sections. For simplicity, we assume the minimum number of massless electric 
quarks, Nfo = 1. We will also assume for simplicity that h = 1; in general, it is not a free 
parameter but is instead determined by various (uncalculable) 0(1) numbers coming from 
Seiberg duality ||. 

The coupling to the messengers is as in ordinary gauge mediation - we couple the field 
$oo (which has both a lowest-component and an F-component VEV) to a vector-like pair 
of messenger fields $ transforming in 5, 5 of SU(5) GUT: 

W mess = A$ 00 ** • (4.1) 

Here A is some dimensionless coupling; since it is actually a non-renormalizable interaction 
in the UV, A will be naturally quite small. We are considering only one messenger pair 
here; the generalization to arbitrary numbers of messengers is trivial and does not affect 
the rough sketch of the physics presented here. 

The overall scale of the model is set by phenomenological considerations. We wish to 
obtain a realistic spectrum of SUSY partners, so we set the gaugino mass to be around 
100 GeV. This is achieved by using the usual formula of ordinary gauge-mediation 

4tt $oo 

Now let us consider the various consistency conditions on the model. In addition to 
the constraints ( |3.20| ), ( |3.21| ) on the SUSY breaking sector, we also have the following 



constraint coming from the requirement of non-tachyonic messengers: 

(A$ 00 ) 2 > AF $00 . (4.3) 

Substituting F$ 0Q = e^oo, this becomes 

A$oo > e/j . (4.4) 

If A were an 0(1) coupling then this constraint would be trivial to satisfy, since $oo 3> A*- 
However, since A comes from a non-renormalizable interaction in the UV, we find that 
A 1 in such a way that (^4.4|) becomes quite restrictive, similar to |35[ . 

To see this in more detail, let the $qo operator be suppressed by a mass scale Mi in 
the electric theory, and the $oo^$ be suppressed by a (possibly different) mass scale M^. 
Then we have 

A A " 
M 2 Mi/i 
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and ( fOD is equivalent to (after substituting $ o ~ "~ with C the combination of color 
and flavor factors and the log appearing in ( |3.6| )) 

M2 <f?«l- (4-6) 



CMi Ae 4tt 

Thus for C ~ (9(1 — 10) there must be a sizeable hierarchy between Mi and M±. In practice 
this turns out to be at least a factor of M2/M1 ~ 10 -4 . This fine-tuning in Mi vs. M\ is 
perhaps the least attractive feature of this model.@ 

We can further manipulate the constraints to derive the allowed ranges of the param- 
eters. We start with (|4.2|) and (|4.5|) , which together imply 

A 2 

efl ~ — ~ 10 5 GeV . (4.7) 



Translating (|3.21|) to the new parameters leads to 



A Mi/i ~ (4tt) 2 ' 1 



where we have taken ft ~ 1 for the upper inequality. Combining this with (177) we obtain 

11 > W 7 GeV , 

A » To4v * 109 GeV ' («) 

Mi > 1= ^ > 10 13 GeV . 

10 15 GeV 3 ~ 

In fact, if we do not want Mi to exceed the Planck scale, then according to (|4.9[) , /i 
is actually bounded above by ~ 10 8 GeV. So we see the parameters can only take values 
in the ranges 

H ~ 10 7 -10 8 GeV , A ~ 10 9 -10 12 GeV , Mi ~ 10 13 -10 19 GeV . (4.10) 



7 It is interesting to note that, if instead of a quartic deformation, the leading non- 

renormalizable term is an even higher-dimension operator, 5W ~ — 2k-3 ^oo? (and assuming the 

M i 

messengers are still coupled through (|4.l|) ) the constraint of non-tachyonic messengers becomes 

trivial for k > 3. For instance, if k = 3 the messengers are non-tachyonic as long as 

Mi A 2 
W 2 > Mf ' 

This is a rather trivial inequality which does not lead to any hierarchy between Mi and Mi. 
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Apparently, the parameters of the model are quite constrained by all the consistency 
conditions, and the result is a relatively predictive version of ordinary gauge mediation! 

Let us conclude this subsection by just mentioning one particularly intriguing choice 
of parameters that satisfies all the constraints (and in fact minimizes the hierarchy between 
M 2 and Mi): 

li ~ 10 8 GeV , A ~ 10 12 GeV, M 1 ~ M pl ~ 10 19 GeV . (4.11) 

This satisfies all the requirements as long as M 2 < 10 15 GeV. Interestingly, this is tan- 
talizingly close to the GUT scale, so perhaps one could imagine a UV completion where 
GUT-scale fields connect $00 to the messengers, but only Planck-suppressed interactions 
generate $g . 

4-2. Phenomenological features of the model 

Now let us briefly mention some phenomenological aspects of our ordinary gauge 
mediation model. For concreteness we will focus on the point (|4.11 ) with M 2 ~ 10 15 GeV. 



However, we would like to emphasize that the statements we will make in this subsection 
are for the most part predictions of our model, not simply consequences of the particular 
choice of parameters. 

In our model, the parameters of ordinary gauge mediation are 

M mess ~ A$ 00 ~ 10 5 -10 6 GeV , 



F ~ \F* 00 ~ 10 10 -10 n GeV 2 

m 3/2 ill ~ 1-10 MeV . 

1 M p i 



(4.12) 



Here we have used A ~ 10 -3 and e ~ 10 -3 . Note that M mess and F come from the same 
source, the R-symmetry breaking scale $00 1 an d therefore there are no extra CP phases in 
the gaugino masses even if doublet/triplet splitting of the messengers is present. Moreover, 
having one and the same source for M mess and F, as we do, guarantees messenger-parity, 
which is necessary for a consistent spectrum of soft terms. 

Since R-symmetry is spontaneously broken, the R-axion is massless in gauge theory, 
up to small non-perturbative contributions. On top of the field theory effects, there is also 



a supergravity contribution to the mass m a . Following j|2| we find that this contribution 
is 




m a ~ ~ 10 GeV . (4.13) 
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This is already well above the astrophysical bound, m a > 10 MeV. 

Note that the mass of the gravitino in (|4.12| ) is always determined by the largest F- 



term in the SUSY-breaking sector. In this model, this is not -F$ 00 but is instead F^ lx ~ 
fi 2 ~ 10 16 GeV 2 . So what we have is a model with light messengers but a heavy (~ 1 MeV) 
gravitino. This combination of features is generally rare in models of gauge-mediation; here 
it is achieved because the hidden sector has a subsector with much larger SUSY-breaking 
than the SUSY-breaking felt by the messengers.! 

A heavy gravitino can have observable consequences at colliders. In general, a grav- 
itino mass of (|4.12| ) implies that the NLSP is sufficiently long lived that it will escape the 
detector. (For a good review of this subject and more references, we refer the reader to @.) 
Thus the standard 77+MET signal of gauge mediation will not be observed. On the other 
hand, if the NLSP is a very long-lived slepton (as would be the case, for instance, if the 
messenger scale is very low or there are N m ~ 5 messengers) then the collider signatures 
of the model could be quite spectacular. 

4-3. Naturalizing the UV model 

To build a natural model (starting from some high scale) with the dynamics described 
above we clearly need to forbid all the super-renormalizable terms in the UV Lagrangian, 
while still allowing for the quartic interactions between the massless electric quarks. In 
addition, as explained in section 3, some quadratic operators must be forbidden for our 
analysis to be reliable. We would also like to forbid the most general couplings between the 
electric quarks and the messengers while allowing for ( |4.1| ); otherwise the messenger SUSY 
masses and F-terms will no longer necessarily come from the same source. Finally, we 
would like to generate the mass scale \i dynamically. As we will see in this section, all four 
of these conditions can be met by "retrofitting the model" [[J6| with an auxiliary gauge 
group and imposing its associated non-anomalous discrete R-symmetry. This approach 
was implemented recently in a similar MGM-type context by ||27[ . 



The main idea of retrofitting is as follows ||36|| . Consider a pure super Yang-Mills 



(SYM) theory SU(N) with the strong scale being A. At energies below A there is a 



The gravitino mass also sets the scale of all the gravity mediated contributions to soft terms. 
These contributions are not degenerate and, in general, violate flavor in the maximal possible way. 
For the gravitino masses considered here, there is no contradiction with precision measurements 
and these contributions are allowed to be arbitrary. For a recent work analyzing flavor aspects of 



combined gauge and gravity mediation see [42]. 
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gaugino condensate which can be used to source masses for some of the electric quarks. 
Furthermore, the SYM theory has a non-anomalous Z 2 ~ R-symmetry, and this can be 
promoted to a gauged discrete symmetry of the entire theory. 

Now let us turn to a concrete realization of the UV theory along these lines. Consider 
the following gauge symmetry: 

SU{3) x SU(N C ) x G SM x Z 6 . (4.14) 

The SU (3) group is a pure SYM theory, Gsm is the SM (or GUT) gauge group and 
SU(N C ) is the gauge group of SQCD. Let Qi denote the Nfi massive electric quarks of 
SQCD and Qo be the only electric quark which is massless in the UV. is generated by 
g such that g 6 = 1. We assign the following charges under "Lq: 

d 2 6 : g~ 2 ; W a W a , Q , Q , * : g 2 ; Qi, Qi : g 3 . (4.15) 

Here W a is the chiral field strength of the SU(3) SYM. Note that these charge assignments 
forbid all the renormalizable couplings. Now let us write down the first few terms in the 
superpotential of this theory: 

W = -^(QoQo) 2 + -^QoQo^ + -^TrW a W a Q 1 Q 1 + . . . . (4.16) 

12 3 

Here . . . stand for operators of dimension 6 and higher, as well as for quartic messenger 
self-interactions, both of which are irrelevant for our discussion. We see that the dangerous 
terms QiQi^^ are forbidden by the discrete R-symmetry. We also see that all but the 
<]>o -type quartic interactions are forbidden.! 

To estimate A we take for simplicity M 3 ~ M p \. For /i 2 ~ 10 16 GeV 2 ( pl| ), we find 



A ~ 10 14 GeV. Hence A ^> A, which is necessary for the consistency of our analysis. 
Gaugino condensation breaks to Z2, which is an unbroken symmetry of the theory. 
Indeed, $00 ~ QoQo, which obtains a VEV, is R-even. Note that this Z 2 R-symmetry is 
naturally extended to the R-parity of the MSSM. 

Lastly, let us briefly discuss the issue of the discrete R-symmetry anomalies. The 
classification and analysis of the various possible anomalies has been discussed, for instance, 



9 We wish to emphasize again that only some of the additional quadratic operators cause 
significant shifts of the vacuum we analyzed (and even in this case it could very well be that 
our qualitative conclusions remain intact), so the fact all of them are forbidden by the symme- 
tries ( |4.14 ) is an added (but by no means necessary) benefit. 
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Ill 



3^,H31H5|. In particular, MB argued that the cubic anomaly does not impose any 



constraints on the low energy effective action. So, we remain with the constraints of 
gravitational anomaly and mixed gauge-discrete anomalies. The mixed anomaly with the 
strong 577(3) is automatically satisfied and the mixed anomaly with SU(N C ) gives 

2N fl + N f0 = 3r , (4.17) 

where r is some integer. 

On the other hand, since we have not specified the R-charges of the SSM fields and have 
not even chosen the matter content ,0 we are unable to calculate the gravitational anomaly 
and the mixed anomaly with the SSM gauge group. The specific charge assignment here 
should be merely considered as an example of how to implement the above mentioned idea 
of naturalization, and we leave a more detailed analysis to the future. 
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Appendix A. Calculating the Runaway Potential at Large Fields 

For the sake of completeness, we will describe in this appendix the RGE evolution 
of the wavefunction of $n for large $oo- A more general derivation of the leading-log 



potential for pseudo-moduli at large fields can be found in f5B|] . 

We will assume that ($00) = Z ^-N f0 f° r simplicity; the generalization to arbitrary 
non-degenerate $00 will be obvious at the end. The wavefunction factors in the Kahler 
potential are governed by the RGEs: 



log^ 11 =2^ A ^ 7<J>11 (g';z), (A.1) 



where 7$^ is its anomalous dimension, and Q is the RG scale. In this model, the relevant 
one-loop anomalous dimensions above the scale hz are 

Nh 2 N fl h 2 N f0 h 2 

7$11 ig^z^z^ ' ^ V ig^z^zi + m^z^z^z^ ' 1 • > 

and below the scale hz they are 

7, - Nh * 7 - N n h " (A3) 

7 *" w^z^z^ ' 7 ^ V i6^z 9ll z^ ■ [A - 6) 

In particular, 7$^ is continuous at the scale hz. Therefore we have 

dlogZ^ f A dQ' gWgj^j o/' A ^ , d{\ogZ^ 1 +2\ogZ Vl ) 

~ 2 ~7V ^ZTTZ =~ 2 nT^u xTZTTZ ■ ( A - 4 ) 



dioghz Jq Q' d log hz Jq Q' 1 9 log hz 

At leading-log order, only the second term in ( A.4| ) matters, and only through the discon- 
tinuity in 7^: 

WF _ [ A d Q' ( N nh 2 , iVfo^ \ , /*" dQ' ( Nfih? \ 

(A.5) 

so 

= , _ e ( fa -«) + ..., (a.6) 

dlog/iz S-k^Z^Z^Z^ 
where ... do not contribute in the leading log approximation. Substituting this into ([A.4Q 
and integrating twice, we conclude that 

Nh 2 Nfoh?f hz\\ 
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where we have dropped the wavefunction factors since they will only contribute at higher 
loop order. 

Now, it is trivial to generalize to arbitrary $00 • For instance, if $00 is diagonal but 
has different eigenvalues, the same arguments as above lead to 



i=l 



Since $00 can always be diagonalized by an SU(Nf )L x SU(Nf )n transformation, the 
full generalization of ( |A.7|) , (|A.8|) must be 



logZ ^ = i6^ Tr y og —q^~ 1 + • • • ' (A - 9) 

which is the result quoted in the text upon substituting Q = h\i. 
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